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Abstract 



Based on Yuen's solutions with radially symmetry of the pressureless density-dependent 
Navier-Stokes in R N , the corresponding ones with elliptic symmetry are constructed by the 
separation method. In detail, we successfully reduce the pressureless Navier-Stokes equations 
with density-dependent viscosity into 1 + N differential functional equations. In particular for 
fti > and K2 = 0, the velocity is built by the new Emden dynamical system with force-force 
interaction: 



with arbitrary constants £, ato and an. We can show some blowup phenomena or global 
existences for the obtained solutions. Based on the complication of the deduced Emden dy- 
namical systems, the author conjectures there exist limit cycles or chaos for this kind of flows. 
Numerical simulation or mathematical proofs for the Emden dynamical systems are expected 
in the future. 
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where the density p = p(t,x) and velocity u — u(t,x) — (ui, U2, ujv) £ R with x = 
(xi,x%, ...,xjv) £ R N with constants 9, Ki, H2 > 0. When 6 — 0, the system is the pressure- 
less Navier-Stokes equation without density-dependent viscosity: 




(1) 



1 Introduction 



The pressureless density-dependent Navier-Stokes equations are 



Pt+V ■ {pu) =0 

p [§- t Ui + u ■ Vu,] =Ki^§r (p 9 V • u) + k 2 V • (p e Vui) for i = 1, 2, N 



(2) 




Pt+V ■ (pu) = 
Ki^f- (V • u) + K 2 Aui for i = 1,2, N 



(3) 
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In fluid dynamics, the Navier-Stokes equations © are the very important mathematical models 
S| and [5]. 

There are intensive papers concerning with the mathematical existence of the Navier-Stokes 
equations with density-dependent viscosity in the recent literature [I], [2], [3], [6], [7] and [l5| . 

In mathematical physics, constructing particular solutions for the Navier-Stokes equations (|2]) 
are an important area for explaining its nonlinear appearance. 

In 2008, Yuen first constructed a class of self-similar solutions by the separation method, for 
the pressureless density-dependent Navier-Stokes equations with radial symmetry [TT]. And the 
separation method is also applied to the other related systems with density-dependent viscosity 
and radial symmetry [T2], [5] and [5]. 

Very recently, Yuen construct a class of self-similar solutions with elliptic symmetry for the 
Euler and Navier-Stokes equations with pressure: 



r pt+v ■ ( P u) -o 

{p^Ui + u-Vui] +K£-p^ Kl £-(V ■u) + K2 £-divu l fori =1,2,...., N w 

with constants K > and 7 > 1. 

The explicit expression for the solutions are 



r N 



n 0fc 

k=l 



H = f i (xi + di) for i = 1,2,...., N 



where 



N 



with s = Yl ^TTTF ' arbitrary constants a > 0, dk and £; 



A- ^ 



and the auxiliary functions Oj = aji(t) satisfy the Emden dynamical system: 

tt, for % = 1,2, ....,iV 



[ II a/; 

v. fc— 1 



<Zi(0) = a i0 > 0, d j; (0) = an 



(5) 



ae 2k s for 7 = 1 

^((-^ + 0)^,0) for 7 >l (6) 



(7) 



with arbitrary constants a^o and an . 

In the paper, we continue to construct the corresponding results for the pressureless density- 
dependent Navier-Stokes equations <j2j> with elliptical symmetry. We successfully reduce the pres- 
sureless Navier-Stokes equations with density-dependent viscosity into 1 + N differential functional 
equations in the following two theorems: 



Theorem 1 To the pressureless density- dependent Navier-Stokes equations (0) with K\ > and 
«2 = in R N , there exists a family of solutions: 

0=^ 

I N 

(8) 

u i =f i (x i +d i ) for i = l,2,...., N 



where 



ae 2k i s for 9=1 

JiS)={ maxff + . 0) for 9 , I (9) 
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N 2 

with s = y~] ft yi , arbitrary constants a > 0, and £; 

k—1 

and the auxiliary functions ai — at (t) satisfy the Emden dynamical system: 

am = y* =1 %4 for i = 1, 2, N 

ai (t) (Ua k (t)) (10) 



«i(0) = a l0 > 0, a,(0) = a,i 

wrf/i arbitrary constants fljo and a^. 
In particular, for £ < 0, 

(1) if all an < 0, the solutions |2p Wow up on or before the finite time 

T = min(-Oio/aii : Oii < 0, i = 1,2, ...,JV); (11) 
^ «/ aZ/ a^i > i/ie solutions (0) esisf globally. 
and 

Theorem 2 To i/ie pressureless density-dependent Navier-Stokes equations 0) tw't/i ki = and 
K2 > in J2 , i/iere exists a family of solutions: 



/(«) 



n a fc 

fc=l 



(12) 



^fo + dO /or t = 1,2, ....,7V 



where 

ae ^2 s /or 0=1 

/(S)= ^max((-^.s + a)^,o) /or ^ 1 ^ 

JV 2 

wii/i s = Y] — iJ , arbitrary constants a > 0, d^ and £; 
fe=l au J 

and i/ie auxiliary functions ai — at (t) satisfy the Emden dynamical system: 



aj(0) = aj > 0, di(0) = a, 



with arbitrary constants a.;o and aii. 
In particular, for £ < 

(1) if some an < 0, t/ie solutions \12\) blow up on or before the finite time 

T = min(-a l0 /aii : au < 0, i = 1, 2, JV); (15) 

(2) if all an > the solutions exist globally. 

Remark 3 For £ > 0, the weak C° solutions |2P and U2\) are the breaking waves with finite mass. 
For £ < 0, the masses of the classical C 2 solutions and h!2\) are infinite. 

Remark 4 Theorem [£| completely covers the result for the pressureless Navier-Stokes equations 
with density- dependent viscosity of Yuen's paper \TTj. 



2 Reduction of the Density-Dependent System 



Similar to Yuen's paper [14] for the Euler and Navier-Stokes equations ((4]), we apply the Yeung 
and Yuen' lemma in |13j and jlOj to obtain the explicit expression for the mass equation: 
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Lemma 5 (Lemma 1 in [10j) For the equation of conservation of mass: 

p t + V • (pu) = 0, 

there exist solutions, 

P 



\ »iW ' »2(t) ajv(t) 



i = |§0=i+di) /or i = 1,2, ....,7V 
with an arbitrary C 1 function f > and > and constants di. 

We can reduce the density-dependent Navier-Stokes equations @ into 1 + N differential func- 
tional equations: 

Proof of Theorem [2j By Lemma [T6l our functions ([5]) satisfy the mass equation 
First, we define the self-similar variable: 



n ai(t) 



(16) 



(17) 



N 

E 



(x fe + d fe ) 2 



(18) 



For the i-th momentum equation of the pressureless Navier-Stokes equations with density-dependent 
viscosity, Ki > and ki — ([2|), we obtain 



AT 



\ dui 



k=l 



d f di \ / cti \ 9 f di 

tt: — (£Ci + di) + — (Xi + di) -— — (xi + di 
dt\ai J \ai J dx.: \a l 



(19) 

-*£|A^P (20) 



AT 



PS 



( di 






(a*)V 



(a) 

(x, + di) + — ^ (xj + d^ 



fc=1 n a k 

k=l 



(21) 



A' 



/O < —(x, t + di 
a,. 



d k f(s) 6 



f(s) 



, a k ( n 

' ' ' II 

fc=i 



6-2 



N 

n a fc 



(22) 



- 2k x 6 



f(s) e - 2 f(s) 



N 

n a k 

k =l 



(23) 



(Xi + di) 24 ) P r 



af II a fc 



fe=i 



(24) 



with the iV-dimensional Emden dynamical system with force-force interaction: 



di{t) 



SI 2^ = 1 a , m , 



ai(t) n a fc (t) 



- for i= 1,2, ...,7V 



ai(0) = aio > 0, di(0) = aii 



(25) 
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with arbitrary constants £, a^ and an. 

Then, we ask the first order ordinary differential equation: 

< 2 «i e / „ w J w or p = 0. 
I /(0) = a > 0, p 

The equation (f2l)|) can be solved by 



(26) 



/(*) 



ae 2k i 
max 



for 9 = 1 

H^ + «)^> o) for^l. 



(27) 



Therefore, the functions © are the solutions for the pressureless density-dependent Navier-Stokes 
equations ©■ 

In particular, for £ < 0, 

(1) if all an < 0, we have 



CLi(t) 



yl~,k=l a, 



■ 



(*) 



> for i = 1,2,..., N. 



It is clear to show that the solutions ([8]) blow up on or before the finite time 

T = min(-a i0 /a l i : an < 0, i = 1, 2, ...,N); 

(2) all Oil > 0, similarly to (1), we have the solutions (|8]) exist globally. 
The proof is completed. ■ 

The proof for Theorem [2] is similar. 
Proof of Theorem [2} For the mass equation, it is the same as the proof in Theorem [T] 

For the i-ih momentum equation ([2]) 2, we have for k\ = and ft 2 > 0: 



JV 



dui s—^ dui 

fe=i k 



k 2 V • (/Vu t 



|- ( — (xi + di)) + \ —{x l + di) ) ^- \ —(xi + dk) 
at V a,- /V a,- / dxi \ a. 



n 2 0ai g_x d 
ai 9 dxi 9 



(28) 



(29) 



(30) 
(31) 



at ( 0i ) I ( ai ) 



K20aj g_ 2 d f(s) 

a % 9 dx % * 

n a k 

fe=l 



(32) 



( Xi + di) a t p J a*af _ /(s) e ~ 2 / (s) 



a. 2 ai 



iv 

n a fc 
fc=l 



(x. t + d l )a i p 
a 2 a; ( II a k 

k—l 



with the iV-dimensional Emden dynamical system: 



a,i(t) = 



a?(t) n a k (t) 



rr for i = 1,2,..., N 



aj(0) = a i0 > 0, dj(0) = a a 



(33) 



(34) 



(35) 
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with arbitrary constants £, a^ and an. 
Similarly, we require: 



^ + /(*) e - 2 /(*)=° or p.O. (36) 



/(0) = a > 0, 
That is 

' ae~^^ s for 9 = 1 



/(*) = < /V fffl-U \ p „ , , (37) 



max 



+ for^l 



Thus, we verified that the functions (| 12[) are the solutions for the Navier- Stokes equations ([2]). 
It is clear to see the following properties of the constructed solutions. 
In particular, for £ < 0, 

(1) if some an < 0, the solutions (|12|) blow up on or before the finite time 

T = min(-a i0 /a l i : an < 0, i = 1, 2, iV); (38) 

(2) if all aji > 0, that the solutions (IT2|) exist globally. 
The proof is completed. ■ 



3 Conclusion and Discussion 

In this article, we successfully reduce the density-dependent Navier-Stokes equations ([5]) into the 
1 + N differential functional equations: 

(f(s), ai (t) fori = 1,2, ...JV), (39) 

to obtain the self-similar solutions with elliptic symmetry. 

Here, the velocity u is constructed by the novel Emden dynamical systems with force-force inter- 
action (fTUl) : 

5i(*) - for i = 1, 2, TV 

a,(t) n a fc (t) (40) 



ai(0) = a i0 > 0, ctj(0) = oa 

or (TT?]) 

/ at{t) = w-xi+l... for i = l,2,...,iV 



(*) (41) 
aj(0) = a l0 > 0, aj(0) = a,i- 



Comparing with the Emden dynamical system ([7]): 

1 = / W C u-i . for i = l,2,....,iV 

°*(A tt \) ( 42 ) 

o,(0) = Oio > 0, ctj(O) = a a 

for the compressible Euler and Navier-Stokes equations without density-dependent viscosity 
the corresponding ones for the pressureless density-dependent Navier-Stokes equations ([2]) are much 
complicated. It is very interesting to investigate this kind of Emden dynamical systems (|40j) and 
(|4Tj) firstly by numerical methods, because it is seen to be difficult to determine the all blowup 
sets, blowup times and asymptotic analysis of the solutions by rigorous mathematics. Then, it is 
natural to conjecture that there exist limit cycles or chaos (especially for the finite mass cases with 
£ > 0) for the Emden dynamical systems PU|) and (HIT) with N > 2. Thus, further researches on 
proving or disproving the conjecture are expected to better understand the constructed flows p]) 
and (|12jl for the system with density-dependent viscosity ([2]). 
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